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Abstract 



0\ 



Soliton strings in mode-locked lasers are obtained using a variant of the nonlinear Schrodinger equation, appropriately modified to 



model power (intensity) and energy saturation. This equation goes beyond the well-known master equation often used to model 
these systems. It admits mode-locking and soliton strings in both the constant dispersion and dispersion-managed systems in the 
(net) anomalous and normal regimes; the master equation is contained as a limiting case. Analysis of soliton interactions show 
^ that soliton strings can form when pulses are a certain distance apart relative to their width. Anti-symmtetric bi-soltion states are 
^ ialso obtained. Initial states mode-lock to these states under evolution. In the anomalous regime individual soliton pulses are well 
'"5 approximated by the solutions of the unperturbed nonlinear Schrodinger equation, while in the normal regime the pulses are much 
T— I wider and strongly chirped. 
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O 1. Introduction 
c/j ■ 

O ' Ultra-short pulses in mode-locked (ML) lasers are the topic 
'55 of extensive research due to their wide range of applications 
J>^ranging from communications 1 1], to optical clock technology 
and even to measurements of the fundamental constants of 
I ^I 'tiature jsll. Although ML lasers have been studied for many 
.years ||4l ll|] it is only recently that researchers have begun to 
fully understand and explore their complicated dynamics. Here 
we discuss pulses associated with a model equation which has 
^ both power (intensity) and energy saturation. This equation 
goes beyond the well-known master equation which is 

contained in the limiting case of small power The master equa- 
tion has gain and filtering terms which are saturated by energy; 
there is an additional nonlinear gain/loss term which grows with 
0^ power (intensity). On the other hand, this model critically dif- 
fers by having the loss term saturated by power (intensity). We 
ghow below that this equation admits soliton strings, including 
\^ anti-symmetric bi-soliton states, in both the anomalous and nor- 
mal regimes and these modes can mode-lock from initial states 
under evolution, i.e. we find these soliton states to be attractors. 
These results are consistent with observations from recent ex- 
periments in both the anomalous |6, 7] and normal |8] regimes 
where in the latter case higher-order anti-symmetric solitons (or 
bi-solitons) were observed. 

A simplified equation without gain/loss terms which is some- 
times used to model ultrashort pulses in mode-locked (ML) 
lasers is the "classical" nonlinear Schrodinger (NLS) equation. 
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The NLS equation in the anomalous regime exhibits both fun- 
damental and high-order soliton solutions. While the funda- 
mental solutions (single solitons) are stationary the higher states 
exhibit nontrivial evolution. In the normal regime the NLS 
equation does not have localized decaying solutions. Indeed, 
the NLS equation has only dark solitons, namely states that 
do not decay at infinity but form on a nontrivial background. 
On the other hand to properly describe ML lasers, in addition 
to chromatic dispersion and Kerr nonlinearity (self -phase mod- 
ulation), saturable gain, filtering and intensity discrimination 
should be taken into account. In this regard, the so-called mas- 
ter equation IH Ht] is an important extension of the NLS equa- 
tion modified accordingly to contain gain and filtering saturated 
by energy (the time integral of the pulse power), while there is 
an additional (sign dependent) cubic nonlinear term that takes 
into account additional gain/loss. The master equation has only 
a small parameter regime where mode-locking to stable soliton 
states occurs |9]; it has not been shown to support higher-order 
states. In fact, for certain values of the parameters this equation 
exhibits a range of phenomena including: mode-locking evo- 
lution, pulses which disperse into radiation, and some whose 
amplitude grows rapidly [10]. In the latter case, if the nonlinear 
gain is too high, the linear attenuation terms are unable to pre- 
vent the pulse from blowing up, suggesting the breakdown of 
the master mode-locking model ||9|]. Thus an improved model 
is desirable. 



If the pulse energy is taken to be constant the master equa- 
tion reduces to a cubic nonlinear Ginzburg-Landau (GL) type 
system. Such GL systems in the anomalous regime have 
been found to support steady high-order soliton solutions ll 111 : 
they also contain a wide range of solutions including unstable, 
chaotic and quasi-periodic states and even blow-up can occur. 
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Furthermore, soliton interaction within the framework of such 
GL systems is compHcated. It is found 1 1 1] that both in-phase 
and anti-phase high-order soliton states are either unstable or 
weakly stable. In either case these states are not attractors 112 1- 
Hence they do not correspond to observations of higher order 
soliton states in a ML laser system which has discrete, nearly 
fixed separations f^. Furthermore, these GL equations do not 
exhibit stable pulses in the normal regime [7J. 

The distributive both power (intensity) and energy saturation 
model which we refer to as the power-energy saturation (PES) 
equation is given in nondimensional form as 
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where il/{z, f) is the slowly varying electromagnetic pulse enve- 
lope, E{z} - \i(f\^dt is the pulse energy, P(z, t) - \tjj'^ is 
the instantaneous pulse power and the parameters g, t, I, e, 6 
are all positive, real constants. The dispersion is represented 
by d{z}, which is taken to either be constant {d{z) - do > 
corresponds to the anomalous regime while d(z) - do < rep- 
resents normal) or a dispersion-managed system (see Section 
|5]l. The dimensionless nonlinear coefficient will be taken to be 
unity (n(z) = 1) in the first part of the article and varying be- 
tween zero and one in the dispersion-managed case, as holds 
in many experiments ifisl iT3l . The first term on the right hand 
side represents gain-saturated by energy, the second is spectral 
filtering-saturated by energy and the third is loss saturated by 
power (intensity). The coefficients e, 6 are related to the sat- 
uration energy and power respectively. If we expand the de- 
nominator of the loss term in a Taylor series in the limit of 
small power and keep only the first two terms, we reduce to 
the master-equation! Furthermore, if one expands the power 
saturation term in Eq. ([T) and keep higher order nonlinearities 
(e.g. fifth order nonlinearities) higher nonlinear GL systems 
result. The energy and power saturation terms in Eq. ([T]l are 
essentially the simplest such forms (inversely linear in their sat- 
uration effects). This type of power saturated loss has been used 
in lumped models to describe experimental ML pulses [ fsl Ht] 
and is widely used as a saturable loss model (e.g. semiconduc- 
tor saturable loss mirrors-SeSAMs). In ifislflol the distributive 
model was employed to show that single soliton states model- 
lock over a large region of parameter space. More recently it 
was shown jlTll that dispersion-managed models with power 
(intensity)-energy saturation are in good agreement with exper- 
imental results in mode-locked Ti:sapphire lasers. It is therefore 
important to study the above distributed model. 

Dimensional values associated with a typical Ti:Sapphire 
laser system are [|l8(] : yS" - 60fs^/mm is the group velocity 
dispersion, 7=1 /cmMW the nonlinear coefficient Pt - 5MW 
the characteristic power, E, - P,t, the characteristic energy, 
Zt - 2mm the corresponding characteristic-nonlinear length 
(Zt - l/yP«), f* - lOfs the characteristic time, g» - 20dB 
the dimensional gain, I, - 0.151/ mm the dimensional distribu- 
tive saturable loss parameter, Alj - 100 THZ the frequency 
cutoff, fisat = lOnJ and Psat - 2MW are the saturated en- 



ergy and power, respectively. Hence: do = /3"z*/f^, g = g,Zt, 
T = gtZtl iS.aP'tl, e - E,/Esat and 6 - -Pt/Psat- These dimen- 
sional values lead to phenomena which are consistent with the 
results presented in this article. 



2. Constant dispersion systems 

The effects of energy and power saturation are crucial in both 
the anomalous and normal regimes. From ML laser behavior, 
the gain and filtering mechanisms are related to the energy of 
the pulse, while the loss is related to the power (intensity) of 
the pulse. Indeed, passive mode-locking generally utilizes sat- 
urable absorbers which are modeled here by a distributive loss 
term saturated with power. 

In the constant anomalous regime id(z) - do > 0) the sat- 
uration terms help localize an intense pulse preventing it from 
reaching a singular state; i.e. "infinite" energy or a blow-up 
in amplitude does not result. Indeed, if blow-up were to oc- 
cur that would mean that both the amplitude and the energy 
of the pulse are large, hence, the perturbing effects are very 
small thus reducing Eq. ([T]i to the unperturbed NLS equation, 
which admits a stable, finite solution. In addition, when local- 
ization occurs, the perturbing effects are small essentially re- 
ducing the equation to the unperturbed NLS equation. Indeed, 
it has been shown 1 16] that when mode-locking occurs the re- 
sulting individual pulses are essentially solitons of the unper- 
turbed NLS equation; the perturbing influence of these terms 
yields the mode-locking mechanism. 

In the constant normal regime {d(z) = do < 0), three regimes 
are observed [19]: (a) when the loss is much greater than the 
gain the pulse decays to zero, (b) when the loss is again the 
prominent effect but sufficient gain exists in the system to sus- 
tain a very slowly decaying evolution resulting in a "quasi- 
soliton" state and (c) the soliton regime above a certain value of 
gain. This soliton differs from its anomalous dispersion coun- 
terpart in that it does not obey the unperturbed NLS equation, it 
is much wider and strongly chirped. 

For example in Fig. [T] a sohton mode is depicted which 
evolves from a unit gaussian state in the anomalous (do - 1) 
and normal (do = -1) regimes (here t - I - 0.1, e - 6 - I and 
g = 0.5 in the anomalous case, g - 1 .5 in the normal case) 

High-order solitons of the anomalous regime {do = 1 > 0) 
are discussed first. Here solitons are obtained when the gain 
is above a certain critical value, g > I, otherwise pulses dis- 
sipate and eventually vanish [16]. As gain becomes stronger 
additional soliton states are possible and 2, 3, 4 or more cou- 
pled pulses are found to be supported. Here we set r = / = 0.1, 
e - 6 - I and vary the gain parameter g. The value of A^/a, 
where and a are the pulse separation and pulse width re- 
spectively, is an important parameter The full width of half 
maximum (FWHM) for pulse width is used, A^, is measured 
between peak values of two neighboring pulses and A<p is the 
phase difference between the peak amplitudes. With sufficient 
gain (g - 0.5) Eq. ([T]i is evolved starting from unit gaussians 
with initial peak separation A^ = 10 and A^/a = 8.5. The evo- 
lution and final state are depicted in Fig. |2] For the final state 
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Figure 1 : (Color online) Mode-locking evolution for a single soliton evolving 
from a unit gaussian initial state in the anomalous (top two figures) and nor- 
mal (bottom two figures) regimes; states at z = 300 are depicted below each 
evolution. For the normal regime the phase-chirp is shown in the inset. 



we find A^/a ^ 10. The resulting individual pulses are sim- 
ilar to the single soliton mode-locking case 11611 . i.e. individ- 
ual pulses are approximately solutions of the unperturbed NLS 
equation, namely hyperbolic secants. The pulses differ from a 
single soliton in that the individual pulse energy is smaller then 
that observed for the single soliton mode-locking case for the 
same choice of g, while the total energy of the two soliton state 
is higher. This is due to the non-locality of energy saturation in 
the gain and filtering terms. 

To investigate the minimum distance, d*, between the soli- 
tons in order for no interactions to occur we evolve equation 
starting with two solitons. If the initial two pulses are suffi- 
ciently far apart then the propagation evolves to a two soliton 
state and the resulting pulses have a constant phase difference. 
If the distance between them is less than a critical value then 
the two pulses interact in a way characterized by the difference 
in phase between the peaks amplitudes: A(f>. When initial con- 
ditions are symmetric (in phase) two pulses are found to merge 
into a single soliton of Eq. ([Til. When the initial conditions are 
anti-symmetric (out of phase by n) then they repel each other 
until their separation is above this critical distance while retain- 
ing the difference in phase, resulting in an effective two pulse 
high-order soliton state. This does not occur in the pure NLS 
equation as shown in the perturbation section below. 

In the constant dispersion case this critical distance is found 
to be A^ - d* ~ 9a (see Fig. |2]i corresponding to soliton initial 
conditions. Interestingly, this is consistent with the experimen- 
tal observations of Ref. To further illustrate, we plot the 
evolution of these cases in Fig. |3] At z = 500 for the repelling 
soUtons A^/a = 8.9. 
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Figure 3: (Color online) Two pulse interaction when < d*. Initial pulses 
(z = 0) in phase: A0 = 0, Af/a » 7 (top) merge while those out of phase by 
A0 = n with Af /ff » 6 (bottom) repel. Here g = 0.5. 



3. Perturbation theory 

To obtain a more fundamental understanding of these numer- 
ical results we use perturbation theory to find evolution equa- 
tions for the six parameters which define the interactions. These 
are: the pulse heights and velocities rjk and Vk respectively for 
k = 1,2, the distance between pulse peaks A^, and the differ- 
ence in phase A0. We take the subscripts 1 and 2 to denote the 
left and right pulse respectively and A0 to be the phase of the 
right pulse minus that of the left pulse. Similarly, A;; = 772 - /71 
and Ay = V2 - V\. Since the individual pulses are well ap- 
proximated by solitons of the unperturbed NLS equation, the 
pulses are considered to be solutions of this equation which 
vary slowly under the perturbing effects of gain, loss, and fil- 
tering. 

Let us denote the NLS equation with perturbation F[il/] as 



1 9 



(2) 



where 



1 + eE 



1 + eE 



il 



1 +6P 



We write the four parameter family of soliton solutions to the 
unperturbed or classical NLS equation {F - 0) as follows 



cr 
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(3a) 
(3b) 



(3c) 



where ji, V, ^0, ctq define height/width, speed, temporal shift 
and phase shift of the soliton respectively. For the classical NLS 
equation without gain, filtering and loss {F - 0), yu, V, ^0, ctq 
are constant. We also note that rj is directly related to yU by 
fi - jf/2. The effect of the right-hand-side being nontrivial is 
approximated by allowing ju, V, ^0, ctq to vary slowly in z. The 
evolution of these parameters is found by satisfying the follow- 
ing integral conservation laws modified by Fli//] and derived 
from Eq. (|2]i 
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(4a) 
(4b) 
(4c) 

(4d) 



where aU integrals are taken over -00 < t < 00. 

Along with the effects of gain, loss, and filtering, the effect 
of tail interaction between the two solitons is also treated as 
a perturbation. Since solitons are widely separated to leading 
order the full solution may be viewed as the superposition of 
two single solitons iff = iffi + il/2', we take t//i to be the soliton on 
the left and i/'2 to be the soliton on the right. When i/'i is locally 
the dominant term i(f2 is treated as a small parameter and we 
expand the nonlinear term about i/'i . If for the moment we omit 
the F[t//] term the evolution of i/'i is well approximated by 



1 2 



-2 |iAi|V2-'A?«A; 



and similarly an equation for when 1^2 dominates is found to 
satisfy 

= G[l/'2,l/'l] 

Notice that if i/'i and i/'2 satisfy this system of coupled PDEs 
then their sum if/ satisfies NLS. Any contributions from inter- 
action in the gain, loss, and filtering terms would be a higher 
order term. For a more thorough treatment of this method 
for analyzing soliton interactions see Ref. |20] and for un- 
perturbed NLS see Ref. 112 111 . Simplifications result if we as- 
sume that Ay - V2 - Vi is small so that from the definition 
A(p =s -yA^ + Act. We also assume A77 = 7/2 - //i is small and 
approximate 771 and 772 as rj. For all variables the bar denotes the 
mean over the two solitons and A the difference of the right mi- 
nus the left. These assumptions are all satisfied in the numerical 
simulations. 

Including the perturbing terms F we have a system of equa- 



4 



tions coupled in their perturbations as follows 



|<Air-Ai =G[<Ai,iA2] + ^'[-Ai] 



ilA2z + -l//2tt + ItAll (^2 = G[t//2, lAl] + F[il/2] 

both of which satisfy Eqs. Substituting the soliton ansatz 
for and tfr2 we find a system of first order differential equa- 
tions for the slowly varying parameters ijk, Vk, ^ok, and crok- 



Subtracting the equations for and we find 



dz 



dl^CTi) 

dz 



which means that A^o, Aery are stationary. Then, from Eqs. Q 
we have A0 x -VA^ + Act and A^ = £{AV)dz + A^o and by 
differentiating we arrive at the following approximate evolution 
equations 



d^ 
dz 



= Ay + A^o,z = Ay 



dA(h 

— - = -y^A^ - yA^, + Ao-^ 

dz 

= -yAy + (^A77 + yAy) + Ao-q,, = fjAri 



The evolution equations are thus found to be 

Sxirik, Vk) + (-1) V sin(A0)e-''^f 
Slim, Vk) + (-1)*^'4773 cos(A0)e-''^f 
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Here a and b are the complex conjugate roots of the polynomial 
+2(1 + 2Tf)x +1=0. The energy for the two NLS soli tons 
is £0 - 477. Note a and b can be chosen to be either root of the 
quadratic equation. A more detailed derivation of these equa- 
tions may be found in f23, Es'], along with an analysis of how 
interactions of solitons in Eq. ^ differ from those in the unper- 
turbed NLS equation. A comparison of the numerical results to 
those found from the asymptotic analysis are seen in Fig. |4]for 
initial conditions 7710 - 7720 = 3.3, A^o = 3 (A^o/a - 6), and 
A^o = 37r/4. 

In Eqs. (Is) the terms and 52 are the contributions from 
gain, loss and filtering, while the other terms come from the 
tail interactions. Considering just the effects of gain, loss and 




Figure 4: (Color online) Numerical simulation (top), asymptotic analysis (mid- 
dle) and comparison of the numerical simulation and the asymptotic result at 
z = 300 (bottom). Here ^> = 0.6. 
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filtering, the velocities have a single stable equilibrium Vi - 
V2 = for any positive value of //i and ri2 (negative values of 
771 or 7/2 are only a phase shift from their positive counter parts), 
so we look at the dynamics of rji and 772 for Vi and V2 at this 
equilibrium. Fig. |5]shows the phase plane for gain values above 
and below the threshold for stable two pulse solutions. In both 
cases there exists a stable equilibrium at 77^ - for A: = 1 or 
2 which amounts to a reduction to the single soliton case and 
an unstable equilibrium when both 771 = 772 = for any g > I. 
Also, in both cases we see an equilibrium at ri\ - ri2 - t]*, 
however, for g < it is found to be unstable and for g > gc, it 
is stable; gc ~ 0.45. 

In Fig. |5] typical cases are shown. Here g = 0.4 (top) is 
unstable and g = 0.6 (bottom) is stable. The 77* found from 
perturbation theory agrees with the height found numerically to 
several decimal places. We also note that there are two other 
two equilibria found for g - 0.6 which are unstable. 



A0() ^ the phase difference will evolve to tt and after some 
initial oscillations the pulses will repel. 
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Figure 6: (Color online) Evolution of the phase for );io = 7/20 = 3.3, Afo = 3.0 
(Afo/o' ~ 6) and several choices of initial phase dilference. Here g = 0.6. 
A<f> = nis shown with a dashed line. 
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Figure 5: (Color online) Typical phase planes for ^ and ^ with Vi = V2 = 



for the unstable two soliton case, 
g = 0.6 (bottom). 



: 0.4 (top) and the stable two soliton case, 



When the soliton interaction terms are taken into consider- 
ation, S I and 5 2 are found to still be the dominant terms for 
sufficiently large A^, hence 7711: rf and — > for k - 1,2. 
However, the small contributions from the interaction terms can 
cause small decaying oscillations for large z- Solving the per- 
turbation theory derived system, we find a stable equilibrium at 
ls.4> - n for any 77 as well as an unstable equilibrium for = 
and A77 = 0. See typical cases in Fig. |6l The peak separa- 
tion is further characterized by the phase difference; the pulses 
attract each other when -njl < A(p < njl and repel when 
nl2 < Acf) < 3.7r/2. This is seen from inspecting the sign of 
AVf - -St?-* cos(A0) exp(^A^) for two solitons with zero veloc- 
ity. Thus, for initial phase difference A0o = the pulses will 
eventually collide and combine into a single pulse, while for 



There is not a non-trivial equilibrium for the full Eqs. ([5), 
so the "effective high order soliton state found numerically is 
not a "true" bound state. We will show, however, that A^(z) - 
0(log(z)) so once the solitons are a certain distance d* apart 
they will be moving too slowly to see or measure. By compar- 
ison, for the classical NLS equation A^(z) - 0{z) |21]. This 
difference is illustrated in Fig. |2]for A0o - ^/^- These "ef- 
fective" bound states emerge naturally from the interaction of 
pulses for most initial conditions making them attractors of Eq. 

o. 
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Figure 7: (Color online) Comparison of NLS to PES; initially the solitons have 
Alp = 71 1 A. For Eq. (T) the resulting pulses have phase difference A(j> = n. Here 
g = 0.6. 

Using the fact that 771 = 772 — » 77* and A(p ^ n as z ^ 00 
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we can derive the asymptotic behavior of A^. First we take the 
difference of ^ and ^ and the derivative of ^ giving us 



_ ^ 4 
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t/^A^ t/AV 



dz 



which may now be combined to form a second order differential 
equation for A^ 
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Since the evolution of A^ is seen from numerical simulations 
to be evolving slowly in z we expect that jJtA^ « ^^'^ 
so ^A^ can be neglected as being a higher order term. This 
leaves a first order equation which we solve to get the long term 
behavior of A^ to be 



-7 log 

77* 



677*(1 +4?7*) 



In fact it is found numerically that the divergence is consid- 
erably less than logz. Hence we call these higher states "effec- 
tively bound" because in practice the laser propagation distance 
is bounded. 

As mentioned above, by increasing the gain parameter addi- 
tional pulses may be supported. In this way we can find three 
and four soliton states that emerge from unit gaussian initial 
states as shown in Fig. |8] In all cases pulses of height 77* 
were taken at an initial separation of A^ = 9a and evolved to 
z - 1000 with less then 0.1 change in peak positions. 
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Figure 8: (Color online) Multiple high-order bound state soliton solutions of 
the anomalously dispersive PES model. 

The depicted solitons are all in phase, however such higher 
order states can be found for any choice of (initial) pair wise 
phase differences. Interestingly, for no interactions to occur the 
critical separation distance between the solitons is again found 
to be d* ~ 9a for the three and four soliton states. For A^ < 



d* the interactions of three and four solitons are more varied, 
but essentially result in either evolving to a lower state through 
the merging of pulses or pulses repelling each other beyond the 
critical distance d*. This behavior agrees with what was found 
for the two sohton state. 

4. Normal Higher Order Solitons 

Next we briefly mention some results for the constant normal 
dispersive case: c/q = -1 < 0. As mentioned above, individual 
pulses of Eq. ([TJ in the normal regime exhibit strong chirp and 
cannot be identified as the solutions of the "unperturbed" NLS 
equation. Indeed, the classical NLS equation does not exhibit 
decaying "bright" soliton solutions in the normal regime. If 
we begin with an initial gaussian, iff{Q, t) - exp(-f^), the evo- 
lution mode-locks into a fundamental soliton state; see Fig. [T| 
bottom two figures. These figures clearly exhibit the mode- 
locking evolution and the significant chirp of the pulse. On the 
other hand we can obtain a higher-order anti-symmetric soliton, 
i.e. an anti-symmetric bi-soliton, one which has its peaks am- 
plitudes differing in phase by n. Such a state can be obtained 
if we start with an initial state of the form if/{Q, t) - f exp(-f^) 
(i.e. a Gauss-Hermite polynomial). The evolution results in an 
anti-symmetric bi-soliton and is shown in Fig. |9] along with a 
comparison to the profile of the single soliton. This is a true 
bound state. Furthermore, the results of our study finding anti- 
symmetric solitons in the normal regime are consistent with re- 
cent experimental observations LSJ. 




Figure 9: (Color online) Evolution of the anti-symmetiic sohton (top) and the 
anti-symmetric (bi-soliton) superimposed with the relative single soliton (bot- 
tom) at z = 1000. Here g = 1.5. 

It is interesting that the normal regime also exhibits higher 
soliton states when two general initial pulses (e.g. unit gaus- 
sians) are taken sufficiently far apart. The resulting pulses. 
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shown in Fig. \TU[ individually have a similar shape to the single 
soliton of the the normal regime with lower individual energies. 
These pulses if initially far enough apart can have independent 
chirps which may be out of phase by an arbitrary constant. We 
again find that d* w 9q' is a good estimate for the required pulse 
separation just as in the constant anomalous dispersive case! 
These pulses are "effective bound states" in that after a long 
distance they separate very slowly (too slowly to measure). 
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Figure 10: (Color online) Symmetric two soliton state of the normal regime 
with in phase pulses. The phase structure is depicted in the inset. Here g = 1.5. 

Additional higher order 3,4,... soliton states can also be 
found. We will not go into further details here. 



5. Dispersion managed systems 

Standard dispersion managed (DM) solitons, like their con- 
stant dispersive counterparts, are obtained from the PES model 
over a wide range of anomalous dispersion ifisll . This is im- 
portant since recent mode-locked laser experiments have been 
conducted in dispersion-managed regimes ifisi Isl [itIi . In cor- 
respondence with Ti:Sapphire DM laser systems we allow the 
dispersion to vary in z as well as introducing nonlinear manage- 
ment in the form of the function n{z) multiplying the nonlinear 
term. Equation ([T]), is used but we change the symbol of the 
envelope from i// - i//(z, t) to u = u{z, t) to distinguish the two 
cases (constant dispersion and DM, respectively). The effect 
of dispersion management is obtained by splitting the disper- 
sion d{z) into two components d{_z) - do + A(z/Za)/Za where 
< Za « 1 is the dispersion-map period. Hence d(z) is large 
and periodic. If d(z) was only (9(1) then the multi-scale av- 
eraging method employed below would result in the constant 
dispersive case discussed earlier. The path averaged disper- 
sion is do and A(z/Za) is a rapidly varying function with zero 
average which we define as follows: A(^) = {-Ai,0 < ^ < 
1/2; A2, 1/2 < ^ < 1}. Here we consider the case of positive av- 
erage dispersion (do > 0). We define the map strength ,s- = Ai /2 
which gives a measure of the variability of dispersion around 
the average. We take do = l,Za = 0.1 and vary s. The effect of 
the nonlinear management is to turn the nonlinearity "off and 
on", i.e. n(0 = {0,0 < f < 1/2; 1, 1/2 < f < 1}. For exam- 
ple, in a Ti:Sapphire laser the nonlinearity is negligible in the 
anomalous regime where one has a prism pair that compensates 
for the normal dispersion in the crystal. The averaged equation 



is derived using the method of multiple scales and perturbation 
theory [24]. The variation in dispersion occurs on the short dis- 
tance scale ^ = z/Za and the pulse envelope evolves on the long 
scale Z - z- 

The method proceeds by expanding u in powers of Za- 

u{(, Z, f) = M<"'(f , Z, f) + Zau'-^\(, Z, f) + 0(zl) 

and substituting this into Eq. ([TJ to obtain a series of equations 
by relating terms by powers of Za- At O(Za^) 

5m<"> A(0 52m(0) 

which can be solved using Fourier transforms to arrive at 
m'°^(^,Z,w) = f/o(Z,w)exp 



where C(^) = Ai^')d^', Uo{Z, oj) = m<"'(^ = 0, Z, w), and the 
Fourier transform pair is defined as 



m(w) - T{u{f)] 



4t)e"^'dt 



00 

1 r^°° 

u{t) = r'^{u(Lo)] = — u{Lj)e-"^'daJ 

Thus Mo separates into a slowly evolving envelope Uo and 
fast oscillations due to changes in the local dispersion. The 
equation for Uo is obtained by imposing secularity conditions 
on the (9(1) terms. 



.dUo do 2 A 
I— (jJ Uo + 

dz 2 



Jo 



exp 

r{F[M(°']})flf^ 



0, 



(6) 



This is the averaged, or mean-field equation which we solve to 
find the pulse dynamics. In the case where F - this is known 



as the dispersion managed NLS or DMNLS equation ||24 1 . 

The method of spectral renormalization |25l] can be em- 
ployed to find single mode-locked DM solitons for Eq. Q lllSn . 
Here we initially superimpose two such DM soliton pulses at 
varying peak separations with phase difference A<p - and let 
them evolve to find two soliton "effective" bound states. As a 
criterion for an "effective" bound state, we require the the peak 
separation differ less than O.OSa after evolving 500 units in z- 
Typical examples with comparison of initial vs. final states are 
depicted in Figs. [TT]and[T2] Much like the single soliton case 
the individual pulses are well approximated by the solutions of 
the unperturbed DMNLS equation. We also note that due to the 
nonlocality of the equation, the individual pulses have a smaller 
energy than the single pulse for the same parameters, as was in 
the case of constant dispersion. 

As is indicated in the figure, the minimum initial distance 
d* X 9a no longer holds for the dispersion (and nonlinearly) 
managed PES equation (DMPES). In the nonlinear managed 
system we find that the critical distance d* ~ la for s - Q, i.e. 
constant dispersion, and d* depends on the map strength s for 
i > 0. Since this change is much more dramatic between s = 



8 



A^/a - 5 




-z = 
-z = 500 



10 




10 15 



Figure 1 1 : (Color online) "Elfective" bound state for two solitons in DM system 
for s = 0. (top) and s = 1.0 (bottom). Here g = 0.6. 



=,°0.4 




"Effective" Bound 
State Distance 



10 



and s - I than s - I and i = 10 we investigate this region more 
thoroughly. In Fig. [T2]the value of A^/a found for the pulses 
to be "effectively" bound are plotted for varying map strengths. 
The general trend is for the needed separation to decrease as s 
increases, however, as can be seen between 5 = 0.1 and s - 1.0 
where more s values were tested this is not a monotonicaUy 
decreasing process. 

For more information on the normal (do < 0) DM case we 
refer the reader to Ref . lfl9ll . 



Figure 12: (Color online) "Effective" bound state for two solitons in DM system 
for i = 10 (top). Numerically found relation between the map strength i and 
the minimum initial distance o / <^ for no interactions to be seen (bottom) after 
z = 500 units. Here g = 0.6. 



6. Anti-Symmetric Bi-Solitons for DM Systems 



We now superimpose two net anomalous DM solitons with 
a TT phase difference. For s < s, 0.25 pulses repel each 
other as was the case for constant anomalous dispersion; we 
note that for the above values of the map strength both of the 
local dispersions being used (do + Ai /za and do + A2/Za) are 
in the anomalous regime. For s > s,, pulses which are taken 
close enough together, i.e. the distance between peak values, 
d < d, ~ 2.5, are found to lock into a bi-soltion state. Examples 
are given in Fig. [13] Pulses taken further apart: d > dt are 
found to repel. These anti-symmetric bi-solitons are the mode- 
locked (due to gain-loss) analog of what was obtained in the 
case of pure DM systems without gain loss 



IE- 



7. Summary 

To conclude, we investigated Eq. dU and found a large 
class of localized solutions including: mode-locked solitons in 
both the constant anomalous and normal regimes, high-order 




Figure 13: (Color online) Bi-Soltion states found for s = 0.3, 0.7 and 5.0. Here 
g = 0.6. The relative phases and the tt jumps at the origin ai'e shown in the 
insets. 
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solitons in the constant anomalous regime and anti-synMnetric 
bi-solitons in the constant normal regime. These results are 
consistent with experimental observations of higher-order soU- 
tons in the anomalous and bi-sohtons in the normal dispersive 
regimes. The dispersion and nonlinear managed system was 
also investigated. Here in the averaged anomalous regime sin- 
gle and higher-order sohton pulses were obtained, including 
anti- symmetric bi-solitons in the net anomalous regime. For 
the constant dispersion case, it is found that when individual 
pulses are initially separated by d* » 9a where a is the width 
of the individual pulse the result is a soliton string. For the DM 
system the results indicate that the high-order solilon strings in 
the DM case can exist in much closer proximity to each other. 
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